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We summarize the extent to which one can use Ward identities (WI) to non-perturbatively improve flavor 
singlet and non-singlet bilinears with three flavors of non-degenerate dynamical Wilson-like fermions. 



It is important to pursue simulations with 
Wilson-like fermions, despite their greater com- 
putational cost than staggered fermions, because 
one avoids the theoretical uncertainty associated 
with taking roots of the fermion determinant. 
One drawback with Wilson fermions is that they 
lead to discretization errors of 0{a), slowing the 
approach to the continuum. This can be amelio- 
rated using an 0(a) improved action, with the rel- 
evant improvement coefficient, csw, determined 
non-perturbatively using the method of Ref. pp. 
One also wants to improve the matrix elements of 
operators, and a technique for doing so has been 
worked out for non-singlet quark bilinears with 
degenerate dynamical quarks pOj. 

Here we sketch the generalization of this im- 
provement program to the realistic case of three 
non-degenerate dynamical quarks, and to all bi- 
linears, including flavor singlets. The inclusion of 
flavor singlets is a necessary part of the technique, 
and is furthermore of phenomenological interest. 

This study is the completion of the work pre- 
sented in [2] , which was incomplete and partly in- 
correct. Many details must necessarily be skipped 
in this short writeup; a complete discussion, in- 
cluding the generalization to two and four flavors, 
will be presented in 

We assume that the following steps in the im- 
provement program have already been taken: (1) 
Csw and ca have been determined by enforcing 
the PCAC relation for non-singlet axial currents 
with degenerate quarks P; (2) Kc is known from 



the vanishing of the "Ward identity mass", so 
that the bare quark masses arrij = l/2Kj — 1/2kc 
can be defined; (3) Non-singlet vector and tensor 
bilinears have been improved by the addition of 
the standard cy and ct terms which are then 
determined by enforcing axial WI in the chiral 
limit |4I5I6| — so that all improved non-singlet bi- 
linears O^-''^)'^ are known in the chiral limit;^ and 
(4) bg has been determined, using the methods 
of Refs. |1I6) or that presented below, and the 
bare coupling go is adjusted, as the mass matrix 
M = diag(mi, m2, is varied, so that the ef- 
fective coupling g^^l + a bgtiM/Nj) is held fixed. 

As already noted in j^j, the generalization to 
non-degenerate quarks introduces many new im- 
provement coefficients. The improved, renormal- 
ized off-diagonal bilinears (O = S, P, V, A, T), 

d^^''^=Zo[l + abotrM+ abo^{mj + mk)]0^^'''>'' 

require one more coefficient, bo, than in the 
quenched case, while the diagonal non-singlet bi- 
linears require another new coefficient, fo- 



q(33) _Q{kk) _ 



{l + abotTM){0'-^^^'^-0'^''''^'') 



-abo{mjO^"'>-mkO'^''''^) + afoimj-mk)trO 



For O = S there is also mixing with the identity 
operator, missed in |2j, and discussed below. 

Flavor singlet bilinears require new coefficients. 
In the chiral limit one has 
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(tryl^)-^ = trA^ + acAdf^tvP 



^Notation: O is a 3 X 3 flavor matrix, flavor labels i,j,k 
are not summed, and are implicitly taken to be distinct. 
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= tTVfj_ + a cvd^tvTfj^^ 
= trT^y + acT [df^trV^ 
= trP + agpTr{Ff,,Ff, 

3, 



a "es + tri> + agsTT{Fi_,^Ff^^) . 

Note that gluonic operators are needed to im- 
prove trP and tr^ (but not trT, as erroneously 
concluded in |2]), and the need to subtract the 
identity operator from trS*. Away from the chiral 
limit one has 

t?d = Zoro [(l + adotrA/)tr(0)^+adotr(A/e))] , 

with new coefficients tq, do and do- is 
scale dependent since the singlet axial current 
has anomalous dimension, while the other tq are 
scale independent. For O = S we have implic- 
itly subtracted mass dependent mixing with the 
identity operator, as discussed below. 
Finally, the improved quark masses are: 

trAAf = Z,„ [{l+abjrMM)tr\M + ab„,tT{\M^)] 

trM = Z^rm[{l+admtTM)trM+ad,ntiiM'^)] 

where A are Gell-Mann matrices. Zm is scale de- 
pendent, while rm is not. Note that there is no 
independent /,„-like term. 

The coefficient eg would need to be determined 
to an accuracy of a'^ to improve trS*, which ap- 
pears impractical. However, to calculate hadronic 
matrix elements one must subtract disconnected 
contributions anyway, and this completely cancels 
the es term. Similarly, the 65 contribution can be 
canceled when implementing Wl by subtracting 
disconnected contributions. In this way one can 
avoid the problem, except when calculating the 
vacuum expectation value. The same issue arises 
for the diagonal non-singlet scalars away from the 
chiral limit, and has the same resolution. 

The first step in determining the improvement 
coefficients is standard: one enforces the normal- 
ization of vector charges for a selection of hadrons 
for various quark masses. This determines Zy, 
ry, by, by, fv, dy and dy, i.e. all the vector 
coefficients except cy and cy. The latter do not 
contribute to the charges. 

In Ref. |2] it was claimed that one can learn 
about other improvement coefficients by enforc- 
ing the vector transformation properties of bilin- 
ears. This is incorrect. The essential reason is 



that the lattice theory does not break the vector 
symmetry. The relations that were to be used can 
be shown to be automatically satisfied. 

The next step is to relate the coefficients for 
scalar bilinears and the quark mass, enforcing the 
standard continuum convention ZsZm = 1. This 
is unchanged from Ref. j2j — the new es term does 
not contribute. The result is that all other coeffi- 
cients for the scalar can be expressed in terms of 
those for the quark mass, so that the former are 
redundant. In particular, we find gs = bg/{2g'^), 
and the constraint ds ^ bs + 865, which follows 
from the absence of an fm term. 

We next consider "two-point" axial WI, i.e. en- 
force the axial divergence relations on improved 
quantities while varying quark masses. In addi- 
tion to the off-diagonal currents used in llfjj and 
0, we now use flavor diagonal currents: 



^.p(w)(a;)_f?i^p(fefc)(a;)] j(o)^ [a; ^ 0] . 
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The RHS involves singlet parts of M and P, il- 
lustrating how these cannot be avoided. When 
expanded out in terms of bare lattice operators 
and masses, there are several independent opera- 
tors on both sides, as well as quark-connected and 
disconnected contractions, so this relation will be 
complicated to implement. Nevertheless, it is a 
powerful constraint, yielding eight independent 
combinations of coefficients 

The final step is to enforce the axial transfor- 
mation properties of the bilinears — the so-called 
"three-point" axial WI. Relations of the form 



are the simplest to implement as they involve no 
quark-disconnected contractions. They have been 
used extensively in quenched simulations |ll5j . 
and can be used in the chiral limit to determine 
Cy and ct as noted above. The generalization to 
the unquenched case is straightforward. 

To determine many of the new coefficients, 
however, one must use other three-point WI. Con- 
sider first the transformation properties of the 
singlet axial current: 
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Here the four- volume V contains y but not z. 
Contact terms are avoided by setting fhi+rhj = 0. 
The above relation allows the determination of 
the relative size of different terms in tr^, though 
not the overall normalization. We find that one 
can obtain ca and cIa- Replacing tr^l with tiV 
one can similarly obtain cy and dy- 

Further coefficients result from enforcing 

Afe')trfT~~ 1 — 

This determines the relative size of different terms 
in trT (giving ct and ^t) and the ratio of normal- 
izations of LHS and RHS (giving rx and hT — dx)- 
The scalar case is of particular interest, i.e. 

(<5f ^5 tr[S{y)] J^^^\z)) ^2{P^^^\y)J^^^\z)) 

From this it appears that one can determine gs 
and ds from the relative normalizations on the 
LHS, and rs and bp — ds from the ratios of nor- 
malizations of the LHS and RHS. There is, how- 
ever, a subtlety because of the identity operator 
contribution to trS*. Formally, this contribution 
vanishes since it it is invariant under axial trans- 
formations, and in particular 

{S'^fSJ^^'\z)) = Oia^). 

However, the divergent coefficient oc inval- 
idates this argument for the lattice calculation, 
so that one must explicitly subtract the discon- 
nected contribution 

from the LHS to remove the eg term. 

This provides a new method for determining 
bg = 2gQgs. Note that one can carry out this 
determination in the chiral limit, so that one does 
not need to know b 



<g a prion. 
The final WI giving new information is'^ 
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The RHS contains singlet and non-singlet pieces, 
and so gives new information on improvement co- 
efficients, although at the cost of requiring quark- 
disconnected contractions. We find that, in addi- 
tion to constants determined previously, this WI 

^This was not considered in 



allows the determination oi fx- A new subtlety 
here is that the contact terms to be avoided can- 
not be obtained simply using the prescription for 
equation-of-motion operators given in . 

Combining results from all the WI considered 
above, we find that one can determine all the 
renormalization and improvement coefficients in- 
troduced above except (1) the scale-dependent 
constants Zt, ZsZp^ and rA (whose determina- 
tion requires use of a method like NPR, which 
fixes operator normalizations by some conven- 
tion); and (2) the two scale independent quan- 
tities dA and 6t + dr-^ 

Several issues merit further investigation. For 
example, is it possible to determine dA and 
br + dr using other Ward identities? How do 
the results generalize to two and four flavors (see 
Ref. [Sj)? What is lost if one works at fixed 
which is easier in practice? And to what extent 
does the proposal of Ref. (which requires an 
even number of flavors in its simplest form) re- 
move the need to determine the improvement co- 
efficients? 
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*This corrects the conclusion of 2 tliat there were three 
undetermined scale independent quantities. 



